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ABSTRACT 


The principal technical challenge in designing a millimeter accu- 
racy satellite to support two color observations at high alti- 
tudes is to provide high optical cross-section simultaneously 
with minimal pulse spreading. In order to address this issue, we 
provide, in this paper, a brief review of some fundamental prop- 
erties of optical re t r ore f lectors when used in spacecraft target 
arrays, develop a simple model for a spherical geodetic 
satellite, and use the model to determine some basic design cri- 
teria for a new generation of geodetic satellites capable of sup- 
porting millimeter accuracy two color laser ranging. We find that 
increasing the satellite diameter provides: (1) a larger surface 
area for additional cube mounting thereby leading to higher 
cross-sections; and (2) makes the satellite surface a better 
match for the incoming planar phasefront of the laser beam. 
Restricting the retroref lector field of view (e.g. by recessing 
it in its holder) limits the target response to the fraction of 
the satellite surface which best matches the optical phasefront 
thereby controlling the amount of pulse spreading. In surveying 
the arrays carried by existing satellites, we find that European 
STARLETTE and ERS-1 satellites appear to be the best candidates 
for supporting near term two color experiments in space. 
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1 INTRODUCTION 

A companion article in these proceedings . [Degnan, 1992] demon- 
strates the benefits of utilizing and maintaining ultrashort pul— 
sewidths in performing two color ranging measurements to 
satellites. Unfortunately, even if one starts with a very short 
pulsewidth (e.g. 35 psec) , most of the existing satellites will 
broaden the pulse significantly thereby degrading the precision 
of the differential time-of-flight measurement. This is espe- 
cially true for the high altitude satellites, such as LAGEOS and 
ETALON, which are the primary targets for space geodesy 
applications. The principal technical challenge in designing a 
millimeter accuracy satellite to support two color observations 
at high altitudes is to provide high optical cross-section simul- 
taneously with minimal pulse spreading. In order to address this 
issue, we provide, in this paper, a brief review of some 
fundamental properties of optical retroref lectors when used in 
spacecraft target arrays, develop a simple analytical model for 
spherical geodetic satellites, and use the model to determine 
some basic design criteria for a new generation of geodetic sat- 
ellites capable of supporting millimeter accuracy two color laser 
ranging . 


2 RETROREFLECTOR CHARACTERISTICS 

For normally incident light, a single unspoiled retroref lector 
has a peak (on-axis) optical cross-section a cc defined by 


Occ = P^< 


f 4ji A. 


V 


X‘ 


( 2 . 1 ) 


where p is the cube corner reflectivity, A cc -nR 2 cc is the light 

collecting area of the corner cube, and AnA cc /X 2 is the on-axis 
retroref lector gain. For a circular entrance aperture, the far 
field diffraction pattern (FFDP) of the reflected wave is the 
familiar Airy function given by [Born and Wolf, 1975] 


o(x)=o 

where 
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The Airy pattern consists of a main central lobe surrounded by 
weak rings. The angular half-width from the beam center to the 
first null is given by the first nonzero root of the Bessel func- 
tion Ji which, with (2.2b), yields the formula 


null 


1 . 22 - 


'D t 


(2.3) 


where X. is the wavelength and D cc = 2 R cc is the diameter of the 
retroref lector . 
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At arbitrary incidence angle, the area in (2.1) is reduced by the 
factor [Minott, 1974] 


T l( 0 tae)"-( sin ~ 1 n--/2jjitan0 r , / )cose, n 


where Q inc is the incident angle and 0 re/ is the refracted angle 
determined by Snell's law, i.e. 



(2.5) 


where n is the cube index of refraction. The quantity n is given 
by the formula 


H = J l - 2tan‘ 


( 2 . 6 ) 


Thus the peak optical cross-section in the center of the 
reflected lobe falls off as 


° e //( 0 mc) = T l 2 (e in c)o cc 


(2.7) 



INCIDENCE ANGLE, DEGREES 


Figure 1: Normalized peak optical cross-section as a function of 
incidence angle for an unspoiled retroref lector. 
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Figure 1 shows the falloff of optical cross-section with inci- 
dence angle for the two most common retroref lectors - hollow 
(n=l) and quartz (n=1.455). Note that, for a solid quartz cube, 
the optical cross-section falls to half its on-axis value at 
roughly 13° incidence angle and is effectively zero beyond about 
40°. The cross-section for a hollow cube corner falls to half its 
normal incidence value at about 9° and is effectively zero beyond 
about 30°. 

One can further limit the effective incidence angle over which 
the retroref lector responds by recessing the reflector in its 
holder. It can be easily shown that the incidence angle at which 
the retroreflector response is zero is given by 



where d is the depth of the recess. 

3 VELOCITY ABERRATION 

As mentioned previously, the far field diffraction pattern (FFDP) 
of a cube corner with a circular entrance pupil function corre- 
sponds to the familiar Airy pattern consisting of a single main 
lobe surrounded by low intensity rings. If there were no relative 
motion between the satellite and the target, the center of the 
FFDP would fall on the instantaneous line of sight between the 
target and satellite. However, due to the relative velocity 
between the satellite and the target, the coordinates of the FFDP 
are translated. The magnitude of the angular displacement in the 
FFDP is given. by the equation [Minott, 1976]: 

a(/i s ,0 2en ,u)) = a max (/i s )>/ cos 2 a> + r 2 (/i s ,9 2en )sin 2 to (3.1) 

where the maximum value, ct m ax^ i s given by the expression 


a. 



(3.2) 


and 


r(/i,.e 2 .„) = 



/? £ sin8 zen V 

R E + h s J 


(3.3) 


and R e is the Earth radius, g = 9.8 m/sec 2 is the gravitational 
acceleration at the surface, h s is the satellite height above sea 
level, c is the velocity of light, and the angle 


uo = cos _1 [(r xp)-fl] 


(3.4) 
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where f, p, and v are all unity length vectors corresponding to 
the satellite position vector (relative to the Earth center) , the 
line-of-sight vector between station and satellite, and the sat- 
ellite velocity vector respectively. Since Y (h s ,Q zen ) is always 
less than unity, equation (3.1) has an effective "minimum" value 
for a given Q zen , when ca) = ji/2. Thus, 

= a m ax( /l S ) r ( ft S ’ 0 ,.J 

The maximum and minimum angular displacements of the FFDP are 
plotted as a function of satellite height in Figure 2 assuming a 
maximum Q zen of 70°. It should be noted from the figure that the 
angular displacement decreases with altitude and that the maximum 
and minimum values converge for high satellites. At ETALON alti- 
tudes (19,000 Km), for example, the angular displacement is 
roughly constant at about 26 jurad. 



Figure 2 : Maximum and minimum angular displacements in the retro- 
reflector far field diffraction pattern caused by the velocity 
aberration effect as a function of satellite altitude above the 
Earth's surface. 

If the target FFDP is angularly narrow relative to the size of 
the velocity aberration displacement, the receiver will lie on 
the low signal edge of the FFDP or even lie outside the FFDP 
entirely. For example, consider a moving retroref lector whose 
face is normal to the ranging system line-of-sight. The cross- 
section is given by (2.2a) with x = kR cc s\n(a.)~kR cc a where alpha is 
the instantaneous angular displacement caused by velocity 
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aberration. If o t^o-nuiir the reduction in cross section is negli 
gible. However, if a is large, the reflected beam will ,, miss ,, the 
receiver and low or nonexistent signal levels will result. For 
non-normal incidence angles, the reflected FFDP is no longer cir- 
cularly symmetric since the collecting (and transmitting) aper- 
ture of the retroref lector appears as an ellipse to the range . 
receiver. The FFDP peak is again along the instantaneous receiver 
line-of-sight, but the FFDP is now given by the two-dimensional 
Fourier transform of the elliptic entrance aperture of the corner 
cube. The velocity aberration causes the retroref lector response 
to be reduced relative to the peak value given by (2.7). This 
reduction is greater for velocity vectors which are parallel to 
the long axis of the ellipse because of the faster falloff of the 
FFDP with angle in this direction. 

4 RETROREFLECTOR "SPOILING" 

To reduce the effects of velocity aberration in large retrore— 
f lectors, the retroref lector is often "spoiled". The goal of 
"spoiling" is to concentrate more reflected energy into the 
annular region bordered by a max and a min . Ideally, one would like 
to uniformly spread the energy within the annular ring yielding 
an optimum cross-section given by 

O id eal = -C(~) (4-1) 

where the quantity in parentheses is the effective target gain 
and n cc is the solid angle subtended by the annular ring of 
interest. However, conventional spoiling techniques generally 
result in average optical cross-sections which agree only within 
an order of magnitude with the ideal limit described by (4.1). 

Spoiling is usually accomplished by introducing slight variations 
into the cube corner dihedral angles (typically less than two 
arcseconds) . This creates a complicated FFDP which, for an inci- 
dent beam normal to the cube face, breaks the initial single main 
Airy lobe into 2N lobes (where N = 1 to 3 is the number of 
spoiled dihedral angles) distributed within an angular annulus. 
The mean angular radius of the annulus increases linearly with 
the dihedral offset angle from a perfect cube and, from diffrac- 
tion theory [see (2.3)], one expects the effective width of the 
various lobes to depend inversely on cube diameter. 

Each of the 2N lobes originates from a different sector of the 
retroref lector entrance aperture. In fact, the FFDP of each lobe 
is determined by the two-dimensional Fourier transform of the 
projection of the 180°/2N sector of the (assumed) circular retro- 
reflector entrance aperture into a plane perpendicular to the 
line-of-sight between the satellite and the station. The 
distribution of energy within this "annulus" is therefore highly 
nonuniform. Furthermore, the effective area for each lobe is 
reduced to: 
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(4.2) 


/ i cc 

A 'ff- ^ B ^2N 

Substituting the latter expression into (2.1), we obtain an 
approximate expression for the peak optical cross-section at the 
center of one of the 2N lobes for the spoiled retroref lector at 
arbitrary incidence angle 


O pe a*( e . n c-^) = T l 2 (e in c)^^ (4 ' 3 

One can also "spoil" the retroref lector by placing or grinding a 
weak lens onto the entrance face. This approach retains the 
single central lobe of the unspoiled cube corner while reducing 
its peak amplitude and spreading the energy over a wider solid 
angle, D. This yields a peak cross-section given by 

(4 ' 4 

and can be an effective approach when velocity aberrations are 
sufficiently small. 


5 SATELLITE OPTICAL CROSS-SECTION 

As noted previously, the optical cross-section which can be 
achieved with a single retroref lector is limited by the need to 
compensate for velocity abberation effects. Received SLR signals 
can only be enhanced by summing the contributions of several 
retroref lectors. Modern geodetic target satellites (e.g., STAR- 
LETTE, LAGEOS, and ETALON) are all designed to be spherical in 
shape in order to avoid the large pulse spreading caused by 
earlier flat panel arrays when viewed at non-normal incidence. 

The spherical shape also simplifies the modelling of nonconserva- 
tive forces acting on the satellite. 

Satellite array size is largely determined by the satellite alti- 
tude since more retroref lectors are required to achieve reason- 
able signal-to-noise ratios over longer slant ranges. Thus, 
STARLETTE (960 Km), LAGEOS (5900 Km), and ETALON (19,200 Km) have 
diameters of 12, 60, and 129.4 cm and average optical cross- 
sections of .65, 7, and 60 million square meters respectively. 

Let us consider a spherical satellite which is uniformly covered 
with retroref lectors. The density of cube corners, as a function 
of incidence angle, is easily seen to be 

A/(e inc )d0 inc = |sine inc d0 inc (5- 1 ) 

where N is the total number of reflectors on the satellite. To 
obtain a simple expression for the overall target cross-section 
cj, we approximate the sum over all of the retroreflectors within 
the allowed range of incidence angle by the following integral 
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(5.2) 


H 

0 = 0 cc f ctQ inc N(Q inc )r] 2 (Q inc ) 

Jo 

where we have used (2.7). If the retroreflectors are not recessed 
in their holders, ri(0, nc ) is given by (2.4). If their angular 
response is limited by the recess, the variation can be well- 
approximated by the expression 


T1 (0 inc )=l- 


e,„ 


e r 


(5.3) 


where 0 max is given by (2.8). Actually, (5.3) is an excellent 
approximation to (2.4) as well provided we choose 0 max equal to 
.54 rad (31°) for hollow cubes or .75 rad (43°) for solid cubes 
respectively (see Figure 1) • Substituting (5.2) and (5.3) into 
(5.2) and evaluating the resulting integrals yields 



(5.4) 


Let us now examine the validity of (5.4) by substituting LAGEOS 
values. The LAGEOS satellite has a radius R = 29.8 cm and is 
imbedded with 426 retroreflectors (422 fused quartz and 4 germa- 
nium) with a clear aperture diameter D cc of 3.81 cm. Ignoring the 
fact that four cubes are germanium, we choose N = 426 and a value 
of 0 max = .75 rad for solid quartz cubes. We now use a value o cc = 
2.834 x 10 6 in agreement with the input values to the RETRO com- 
puter program as determined during LAGEOS testing and evaluation 
[Fitzmaurice et al, 1977], Substituting the latter values into 
(5.4) yields 

a laceos = 9 . 8 o cc = 2 . 78 x 10 7 m 2 C 5. S ) 


This is roughly equal to the peak value computed by the much more 
detailed RETRO program which showed a range of values between .54 
and 2.7 x 10 7 m^ . Equation (5.4) tends to overestimate the actual 
cross-section because it includes only geometric, and not veloc- 
ity aberration, effects. Equation (5.5) also implies that the 
LAGEOS array cross-section is roughly 9.8 times that of a single 
cube corner at normal incidence. 
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Figure 3: Diagram of a spherical geodetic satellite defining the 
variables used in the analysis. 


6 SATELLITE IMPULSE RESPONSE AND TARGET SPECKLE 

An individual retroreflector responds as a point source and hence 
does not spread the laser pulse in time. However, with a typical 
array of retroref lectors, the laser pulse arrives at the "reflec- 
tion center" of each retroreflector at a slightly different time 
leading to a broadening of the received pulse [Degnan, 1985] . The 
location of the "reflection center" for an individual solid cube 
corner reflector is given by [Fitzmaurice et al, 1977; Arnold, 
1978] 


A R(Q inc ) = nL 



sin 0, 


n 


nL cos0 re/ 


( 6 . 1 ) 


where A/?(0 inc ) is measured from the center of the front face of 
the cube corner to the reflection point, L is the vertex to front 
face dimension, n is the refractive index of the corner cube 
material, 0,- nc is the angle of incidence, and 0 re/ is the corre- 
sponding refraction angle. From Figure 3, it can be seen that a 
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cube at an angle 0 mc to the incident wave produces a time delay, 
relative to the surface of the satellite closest to the ranging 
station (B inc = 0) , given by 

p 2 R,( 

/M<iQ iric ) = ^{R s -[R s -AR(B inc )]cosQ inc } = —(y\-cosQ inc 

where R s is the satellite radius. It should also be noted that 
the differential delay between target reflection points also 
introduces a random phase delay between individual reflectors. 
Thus, if the temporal profiles from multiple cubes overlap at the 
range receiver, the electric fields will interfere . with each 
other in a random way from shot to shot resulting in target 
"speckle". On average, however, the return waveform from the sat- 
ellite should behave as if each of the retroreflectors is an 
incoherent source. This was an implicit assumption in our 
derivation of target optical cross-section in Section 5. 

In the same spirit, the time-averaged satellite impulse response 
can be estimated by summing the weighted (incoherent) returns 
from each of the retroreflectors. Using the simple model for a 
spherical satellite introduced in Section 5, the impulse response 
can be shown to be 


1 - — COS0 


re/ 


( 6 . 2 ) 


Z Jo 

where the geometric weighting factor is given by (2.4) or 
(5.3) ,Ai(0; nc ) is given by (6.2), and the delta function 
6[f - At(Q inc )] represents an infinitely short laser pulse waveform 
incident on the satellite. From (6.2), we see that the delta 
function is nonzero only when the condition 


1 -x 


cos0(t,e,n)-- r , 

1 _e v 1 _ ^ L — • — J 

holds where we have defined the new variables 


(6.4) 





max 


6 = 


nL 

~R~ S 


1 -cos0 max 



(6.5) 


The variable x is a normalized time, expressed in units of the 
roundtrip transit time from the surface of the satellite to the 
center and back, and e is the ratio of the optical depth of the 
cube to the satellite radius. The minimum and maximum values of x 
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are determined by setting 0 equal to zero and 0 max respectively in 
(6.4) and solving for x. The total pulse duration, measured at 
the baseline, is given by A t = t max -e. 

From (5.3) and (6.3), the satellite impulse response can now be 
expressed as a function of the variables x, e, n, and 0 max , i.e. 


/ (x,e,rc,0 max ) 


N . 

o cc — sin0(x,e,n) 


1 - 


0(x , e, n) 

0 maY 


( 6 . 6 ) 


where 0(x,e,n) is defined by (6.4). In the limit of large satel- 
lite diameters (e->0), (6.4) reduces to the simple form 

0(x , 0 , n) = cos" 1 ( 1 - x) (6.7) 

and (6.6) becomes 


/ (x , 0 , n , 9 max ) 



cos' 1 ( 1 - x) 


T2 


( 6 . 8 ) 


The quantity £ is typically small and, for nonzero values of e, 
(6.4) can be easily solved by iteration using (6.7) as a starting 
point, i.e. 


[cos0(x , e , ft)];*! 


1 -x 


1 - eJ 1 


1 ^ |" [COS0(T, £, n)], "J3 

7 2+ l n J 


(6.9) 


until it converges. 

As an illustration, let us use (6.6) to estimate the impulse 
response of the LAGEOS satellite. Substituting n = 1.455 (fused 
silica), L= 1.905 cm, and R s = 29.8 cm into (6.5b), we obtain a 
value e = .093. We recall that, for solid cube corners, we can 
use a value 0 max = .75 rad. Now, using (6.6) and (6.9), we obtain 
the plot of the LAGEOS impulse response shown in Figure 4(a). The 
profile shows the characteristic fast rise and long tail of the 
LAGEOS response. Furthermore, if we compute a center-of-mass cor- 
rection from the centroid of this impulse response profile, we 
obtain a value of 250.2 mm which is in excellent agreement with 
the accepted LAGEOS value of 249 ± 1.7 mm [Fitzmaurice et al, 
1977]. 

The impulse response in the large satellite limit, given by 
(6.8), is shown in figure 4(b) for the case of solid quartz and 
hollow cubes. Note that both the temporal width of the reflected 
pulse and the target optical crossection is smaller for the hol- 
low cubes than for the solid cubes because of their smaller field 
of view. 
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Figure 4: (a) Impulse response of the LAGEOS satellite as com- 

puted by our simple analytical model, (b) Impulse response in the 
large satellite limit for both hollow and solid quartz cubes. 

7 FEASIBILITY OF MILLIMETER ACCURACY SATELLITES 

We will now demonstrate that, in order to achieve high optical 
cross-section simultaneously with minimal pulse spreading with a 
spherical satellite, we must increase the satellite diameter 
and/or retroreflector density and simultaneously restrict the 
response to retroreflectors within a relatively small solid angle 
on the satellite surface about the station line of- sight. 
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The total time duration of the reflected pulse (0% to 0% inten- 
sity points) can be determined from (6.5), i. e. 


2 R s 2 R, 

A< = — (T ma X -G) = — 


n 

~ . 2°max 

2sin — — + g 


L sin 2 0 max , 

cos e maxA / 1 — 2 1 


(7.1) 


which, in the limit of small maximum incidence angles, reduces to 


R s o 

a*~— e 2 max 


1-6 l+-= 


n 


(7.2) 


In the same limit, Eq. (5.4) for the satellite optical cross- 
section reduces to 


o~ 


24 


e: 


(7.3) 


As mentioned earlier, the angular response can be restricted by 
recessing the retroref lectors in their holders. Substituting 
(7.2) into (7.3) yields 


o = 


o cc N 


cAt, 


24 


(7.4) 


We can now express the total number of retroreflectors as 

N = 4(3^| (7.5) 

' D 2 


where p is a "packing density" ( = .435 for LAGEOS) which repre- 
sents the fraction of total surface area occupied by the cube 
faces. Substituting (7.5) into (7.4) yields our final result 


o cc fiRjCAt 


(7.6) 


The product [3/? s cAf in (7.6) quantifies our earlier statement 

that, if we wish to reduce the amount of pulse-spreading A t by 
some factor (via reduction of the retroref lector f ield-of-view) , 
we must increase the retroreflector packing density -satellite 
radius product (|3/? s ) by the same factor to retain a similar tar- 
get cross-section. 


8 CONCLUSION 

As mentioned previously, the principal technical challenge in 
designing a millimeter accuracy satellite to support two color 
observations at high altitudes is to provide high optical cross- 
section simultaneously with minimal pulse spreading. Increasing 
the satellite diameter provides: (1) a larger surface area for 

additional cube mounting thereby leading to higher cross- 
sections; and (2) makes the satellite surface a better match for 
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the incoming planar phasefront of the laser beam as in Figure 3. 
Simultaneously restricting the retroref lector field of vie * 
by recessing it in its holder) limits the target response to the 
fraction of the satellite surface which best matches the optical 
phasefront thereby reducing the amount of pulse spreading. 

For near term experiments, the small radius of STARLETTE makes it 
an attractive target for testing and evaluating two color syste 
or for testing atmospheric models. Furthermore, its low altitud 
(960 Km) and moderate target cross section results in relatively 
high received signal levels. AJISAI, also in a relatively l°w 
1375 Km orbit, consists of small clusters of retroreflectors p 
arated by large reflecting panels and has comparable signal 
strength to STARLETTE. Unfortunately, the satellite is quite 
large and simultaneous returns from several retro clusters 
results in a complicated satellite signature [Prochazka et al, 
19911. From Figure 4, LAGEOS spreading is in excess of 150 pico 
seconds FWHM in agreement with [Fitzmaurice et al, 1977] 
although, with sufficiently short laser pulses (<50 psec) , 
individual retro rings should be resolvable via streak cameras at 
certain satellite orientations. (Note that the simple satellite 
model presented here gives an average response over the full 
range of satellite orientations and shows none of the structure 
expected from a particular orientation) . The LAGEOS pulse spread 
ing combined with relatively low signal returns, measured at the 
few to several tens of photoelectron level for ^st systems, 
would make the necessary differential timing very difficult [Deg- 

nan, 1992] . 

Another useful target for two-color system evaluation is th ® 
recently launched European Earth Remote Sensing satellite, ERS-1. 
It flies at a relatively low altitude (<800 Km) and has a small 
compact target consisting of one nadir-viewing retroref lector 
surrounded by a uniformly spaced ring of eight identical cube 
corners at a nadir angle of 50°. Approximate modelling of this 
satellite by the author indicates that sharp returns consisting 
of one or two peaks (well separated) can be obtained from most 
viewing angles as illustrated in Figure 5. At nadir angles 
between 0° and 15°, the nadir viewing cube is dominant whereas, 
for nadir angles between 30° and 70°, the ring provides a sharp 
return at virtually all azimuthal angles. At nadir angles between 
about 15° and 30°, there is some overlapping of returns and pulse 
distortion. 
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ERS-1 TEMPORAL RESPONSE (50 psec PULSE) 
AZIMUTH ANGLE 



Figure 5: Approximate (geometric) response of the ERS-1 satellite 
to a 50 picosecond pulse as a function of nadir and azimuthal 
angle . 
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